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Ultrathin FeSe films grown on SrTiOs substrates are a recent milestone in atomic material 
engineering due to their important role in understanding unconventional superconductivity in 
Fe-based materials. Using femtosecond time- and angle-resolved photoelectron spectroscopy, 
we study phonon frequencies in ultrathin FeSe/SrTi 03 films grown by molecular beam epitaxy. 

After optical excitation, we observe periodic modulations of the photoelectron spectrum as a 
function of pump-probe delay for 1 unit cell, 3 unit cell, and 60 unit cell thick FeSe films. The 
frequencies of the coherent intensity oscillations increase from 5.00 ± 0.02 to 5.25 ± 0.02 THz 
with increasing film thickness. By comparing with previous works, we attribute this mode to 
the Se Aig phonon. The dominant mechanism for the phonon softening in 1 unit cell thick 
FeSe films is a substrate-induced lattice strain. Our results demonstrate an abrupt phonon 
renormalization due to a lattice mismatch between the ultrathin film and the substrate. 

Ultrathin films regularly exhibit exotic electronic properties. Recent examples include the gi¬ 
ant electron mobility in graphene,^ and the indirect-to-direct bandgap transition in metal chalco- 
genides.l^ Notably, a superconducting gap opening at high temperatures was discovered in one unit 
cell thick FeSe films on SrTiOs substrates (1 UC FeSe/STO).l2J0 The gap opening temperature 
(Tc) of nearly 70 K exceeds the of bulk FeSe by one order of magnitude,!^ and sets a new record 
for Fe-based superconductors. 1 UC FeSe represents the fundamental building block of single 
crystal Fe-based superconductors and can provide critical information for understanding the entire 
material class.^^® 

Extensive studies on Fe-based superconductors have found that the Fe-chalcogen or Fe-pnictogen 
bond angle strongly impacts magneticISHni and superconducting properties.In FeSe/STO sys¬ 
tems, these properties are also strongly dependent on the FeSe film thickness. EJClit 

is thus essential 

to study lattice properties of FeSe/STO as a function of the film thickness. 

Femtosecond time- and angle-resolved photoelectron spectroscopy (trARPES) is a powerful 
technique to investigate how electronic structures respond to lattice dynamics.^^^^ With optical 
excitation, this technique drives coherent phonon oscillations and determines phonon frequencies 
with ~ 0.01 THz precision.E^In Ee-based superconductors, coherent excitation of the generic Aig 
phonon mode directly modulates the Ee-chalcogen or Ee-pnictogen bond angle,which 
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can be related to transient tuning across the complex phase diagram.^ 

In this Letter, we combine a state-of-the-art trARPES setupA^Hl molecular beam epitaxy 
(MBE)l^to study EeSe/STO systems. We observe coherent phonon oscillations in 1 UC, 3 UC, and 
60 UC FeSe films. The phonon frequency evolves from 5.00±0.02 THz for the 1 UC film to 5.25 ± 
0.02 THz for the 60 UC film. By comparing to a Raman scattering experiment,!^ we attribute this 
mode to a Se Aig phonon. Our phonon frequency calculation based on density functional theory 
(DFT) suggests that the Aig phonon softening in 1 UC films is predominantly due to a substrate- 
induced lattice strain. 

Details of the experimental setups for trARPES and MBE are discussed elsewhere .1^!^ We an¬ 
neal 0.05%-wt Nb-doped SrTiOs substrates at 800°C for 15 minutes. Ultrahigh purity selenium 
(99.999%) and iron (99.995%) are then evaporated to the SrTi 03 substrates at 380°C, resulting 
in high quality thin films checked by reflection high-energy electron diffraction. The films are 
transferred to the trARPES system via a vacuum suitcase with a base pressure of 5 x 10“^® Torr. 
For trARPES, we use a Coherent RegA Ti:Sapphire laser system operating at 312 kHz repetition 
rate. We excite the thin films using ~50 fs, 1.5 eV pump pulses, and probe the ensuing dynamics 
by photoemitting electrons at variable time delays using ~100 fs, 6 eV pulses. The energy reso¬ 
lution is ~ 22 meV. For all measurements we use an incident pump fluence of 0.54 mJ.cm“^, and 
maintain an equilibrium sample temperature of 20 K. We do not resolve dependence of the phonon 
frequencies on the pump fluence within the range of 0.07 ~ 0.54 mJ.cm“^. DFT calculations were 
performed using the Quantum-Espresso package using PAW pseudopotentials.!^ For comparison, 
calculations were done using the PBEP^and PBEsol!^ functionals. We used allxllxV momen¬ 
tum grid; the wave function and density cutoffs were 60 Ry and 600 Ry, respectively. 

We show trARPES spectra near the Brillouin zone center for 1 UC, 3 UC, and 60 UC FeSe/STO 
in Fig. [T]^a)~(c), where the photoemission intensity is plotted as a function of electron energy, 
momentum, and pump-probe delay. Before time zero, the energy-momentum cut gives an ARPES 
spectrum in equilibrium. For the 1 UC film, two bands are identified near 80 meV (Band 1) 
and 200 meV (Band 2) below the Fermi level (Fp). We do not observe additional features at 
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(a) 1 UC FeSe/SrTiOg (b) 3 UC FeSe/SrTiOg 


(c) 60 UC FeSe/SrTiOa 



1.0 - 

o 


1.0 - 

o 

\.x-~ - 

1.0- 

o 


^ 0.8- 


^ 0.8- 


^ 0.8- 


0.6 - 

^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

0.6 - 

1—1—1—1—1—1—1—1—1—1—1—1—H 

0.6- 

1—^^^^^—1—1—^—1—1—1—l—^ 

^ 1 - 

I I I 

_ i 1 

^ 1 - 

1 1 1 

A i i i i 1 A A aA 

^ 1 - 
^ n_ 

1 1 1 

o u 
^ -1 - 

^^^^^^^^^^^^ 

o U 

—1—1—1—1— I I—1—1—1—1—H 

o u 
^ -1 - 

vv\r\r\/\A^\A^^'^^ 

<\ 

o 1 

0 - 

< -1 - 

1 1 1 

< 

o 1 

0 

< -1 - 


< 

_p 1 

0 

< -1 - 

1 1 1 1 1 1 1 1 1 1 1 1 

— ^^^^^^^^^^^ — n 


0 1000 2000 0 1000 2000 0 1000 2000 


Delay (fs) 


Delay (fs) 


Delay (fs) 


Figure 1: Extraction of coherent oscillations from the trARPES data on EeSe thin films. (a)~(c) 
3D visualization of the photoemission intensity as a function of energy, momentum and pump- 
probe delay for (a) 1 UC, (b) 3 UC, and (c) 60 UC FeSe thin films on Nb-SrTiOs. We define 
two integration windows for Band 1 (blue) and Band 2 (red). (d)~(f) Photoemission intensity 
dynamics (red and blue dots) corresponding to the two integration windows as defined in (a)~(c). 
Iq stands for the intensity before time zero. We use a lOth-order polynomial to account for the 
smooth backgrounds (solid black lines), and plot the fitting residuals which reveal oscillations 
with periods 200 fs. 
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lower binding energies as reported in the literature.l^ The observed bands are shifted up towards 
£'f when increasing the film thickness from 1 UC to 3 UC, which is consistent with a previous 
ARPES investigation and signifies a change in the charge transfer from the substrate to the film.l^ 
According to the literature,!^ there are two hole-like bands in the spectral region of Band 1 in 
multi-layer FeSe films. In our data, a fine splitting near Ep is indeed observed for 60 UC films 
before time zero, yet it is not observed for 3 UC films possibly due to the finite energy resolution. 

At time zero, the optical excitation creates a non-equilibrium electron distribution. For all 
film thicknesses, the spectral peaks for Band 1 near Fp are almost depleted near time zero, while 
those for Band 2 are partially depleted. During the relaxation process, periodic modulations of the 
spectral intensities can be observed for all three thicknesses. As shown in Fig.[T]^a)~(c), we choose 
energy-momentum windows to be centered at these bands, and integrate the intensities within these 
windows. The corresponding intensity dynamics in panels (d)~(f) clearly display oscillations with 
periods of ~ 200 fs (frequencies of ~ 5 THz). We subtract smooth backgrounds by fitting the 
intensity dynamics with a 10* order polynomial. All fitting residuals exhibit these oscillations. 

We present a more detailed frequency analysis in Fig. To highlight small frequency shifts, 
we use dashed lines to trace the oscillation peaks in Fig.|^a). While the oscillation peaks at early 
delays are temporally aligned, the peaks for the 1 UC and 3 UC films at t > 1000 fs are systemati¬ 
cally delayed. This comparison suggests a softening of the oscillation frequency for thinner films. 
Quantitatively, we fit the residual traces with the functional form of Aoexp {—t/x) cos {2nft + <!>). 
The frequency softens from 5.25 ± 0.02 THz for the 60 UC film to 5.00 ±0.02 THz for the 1 UC 
film, in agreement with the trend in the Fourier transforms of the residual traces (Fig.|^c)). Table[T] 
further demonstrates that this softening applies to both observed bands. 

Table 1: Comparison of phonon frequencies extracted from fitting with a modified cosine function 
(Fig. I^a)) for different bands and film thicknesses 


Thickness (UCs) 

Band 1 (THz) 

Band 2 (THz) 

1 

5.00 ± 0.02 

4.99 ± 0.03 

3 

5.11 ±0.02 

5.10 ±0.03 

60 

5.25 ± 0.02 

5.25 ± 0.02 
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Band 1 


(a) 



Delay (fs) 


(b) (c) 



Figure 2: Analysis of the coherent oscillations for Band 1. (a) Direct comparison of the fitting 
residuals for different film thicknesses. A frequency shift is identified by comparing the corre¬ 
sponding peaks of the oscillations (dashed black lines). These oscillations are fitted by a cosine 
function multiplied by an exponential decay envelope (solid pink lines), (b) Comparison of the 
frequencies extracted from cosine fitting, (c) Comparison of the Fourier transforms of the residual 
curves. The frequencies from cosine fitting are marked by dashed lines. 
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Figure 3: Time-dependent dispersion analysis for 60 UC FeSe/STO. (a) Examples of energy dis¬ 
tribution curves (black) near F and the corresponding fits (red) at 1500 fs. (b) Dispersion of Band 
1 at 1500 fs (solid circles) fitted using a parabolic function (solid curve), (c) Rigid energy shift 
(Frigid, bottom hollow circles) as indicated in panel (b). The oscillatory component (AF, top hollow 
circles) is extracted by removing the smooth background approximated by a 4th-order polynomial. 
AF is then fitted by a cosine function multiplied with an exponential envelope to extract the fre¬ 
quency of 5.23 ± 0.02 THz. 

In addition to their intensities, the binding energies of the bands also exhibit ~ 5 THz oscil¬ 
lations. We demonstrate this using the 60 UC data, as the sharp bands observed therein permit 
a robust dispersion analysis. Since no band splitting is observed near Fp after the pump excita¬ 
tion, we fit the energy distribution curves using two Gaussian peaks multiplied with a Fermi-Dirac 
distribution (Fig. |^a)). In Fig. [^b) we demonstrate the corresponding dispersion for Band I at 
1500 fs. These transient dispersions are then fitted using a parabola with a delay-dependent rigid 
energy shift (Fig. |^c)). This energy shift exhibits oscillations of frequency 5.23 ±0.02 THz in 
agreement with the intensity analysis in Fig. 

To determine whether the ~ 5 THz modes originate from FeSe or STO, we perform a control 
experiment on a Au film of a few nm thickness on a STO substrate. As shown in Fig.|^ the data on 
Au/STO displays no detectable oscillations, which suggests that the ~ 5 THz modes on FeSe/STO 
are primarily from FeSe thin films. 

In accordance with a Raman scattering experiment on bulk FeSe,l231 we assign the ~ 5 THz 
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Figure 4: Comparison of trARPES experiments on 1 UC FeSe/STO and Au/STO. We use a 10th- 
order polynomial to account for the smooth backgrounds and extract the fitting residuals (bottom 
plots). 
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modulations to the Se Aig mode. Aig modes are in general the dominant modes coherently excited 
by ultrafast optical excitations.l^This has been demonstrated in a variety of materials such as Bi ® 
Sb,I^Bi 2 Se 3 ,^and Fe-based superconductors 2|2 1 122 | observation of a dominant Aig mode 

is consistent with this general picture. 

To understand the softening of the Aig mode with decreasing film thickness, it is important to 
consider the strong coupling between the first UC of FeSe and the STO substrate.l^JS One critical 
effect due to the substrate is the FeSe lattice strain.*^*^ The in-plane lattice constant a of FeSe 
increases from 3.77 A in bulk crystals!^ to 3.9 A in 1 UC thin films matching the STO lattice con¬ 
stant.!^ As shown in Fig. this lattice strain causes the Fe-Se-Fe bond angle to increase, which 
results in a lower restoring force for out-of-plane Se displacements. Here we use DFT calculations 
to compare the Aig frequencies for 1 UC FeSe with a freely relaxed lattice constant, and the one 
with a constrained lattice constant. We find that the lattice strain produces a 4% Aig phonon soften¬ 
ing, which agrees well with our experimental value of 5% softening. Additional mechanisms such 
as interlayer interactions! 1 9 | 30|3 1 1 interface dipole field ^ ! ^^!^^ may also need to be incorporated 

for a complete understanding. 

Our experiment demonstrates the high potential of integrating MBE sample growth with trARPES 
measurements. In particular, it allows us to resolve the phonon softening of 0.25 ± 0.03 THz 
(1.0 ±0.1 meV) on multiple bands. This result would be difficult to obtain using traditional 
ARPES, and complements Raman spectroscopy by providing band sensitivity. With the capa¬ 
bility of engineering materials on an atomic scale and studying them in non-equilibrium, this com¬ 
bined approach provides a platform to design and characterize electronic and lattice properties of 
nanoscale materials. One limitation is that the photon energy of 6 eV does not allow us to reach 
the electron pocket at the Brillouin zone boundary where the superconducting gap is located.®^ 
This limitation can be overcome using higher energy probe photons.l^S 

The observation of the Aig mode for the chalcogen or pnictogen atoms in the Fe-based su¬ 
perconductors has important implications.In particular, this oscillation transiently mod¬ 
ulates the Fe-chalcogen or Fe-pnictogen bond angle (Fig. |^, which in turn modifies the effective 
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Figure 5: Comparison of different film geometries impacting the Aig mode frequency. The Fe-Se- 
Fe bond angle is larger in 1 UC films, resulting in a smaller net restoring force. 
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carrier density in the Fe 3(1 orbitals.A number of experiments have shown that the bond angle 
is correlated with for a variety of Fe-based superconductorIt is intriguing that the ob¬ 
served Aig phonon frequency is softened in 1 UC FeSe films, which coincides with a dramatically 
enhanced T^. Understanding the relationship between bond angle, phonon frequency, and electron- 
phonon coupling strength will be an important ingredient to identify the microscopic mechanism 
for high Tc in 1 UC FeSe/STO. 
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